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INTRODUCTION 


\ 

Acoustic  propagation  in  the  Arctic  is  intimately 
related  to  the  structure  of  surface  pack  ice.  In  this  paper, 
reflection  loss  and  transmission  loss  are  calculated  based 
upon  statistical  models  of  ice  roughness.  in — Section  4-r- ~ 

statistical  models  of  ice  ridges  are  introduced  and  tech¬ 
niques  for  estimating  the  roughness  spectrum  of  ridged  sur- 
faces  are  developed.  In— Sect  ion- -%$'  a  model  for  estimating 

the  mean  reflection  loss  is  introduced.  The  reflection  loss 
formulas  are  based  on  the  statistical  models  for  roughness^ 

Cl 

introduced  in  Section  IV  Ivy.  Section-  Jv  transmission  loss 

P;»  r*  ,jr->  i  c 

calculations  using  High-Angle (ARCHAPE)  in  a  simulated 
Arctic  environment  are  compared  with  data.  The  results 
demonstrate  that  rough  surface  scattering  from  random-depth 
ice  keel  structures  can  account  for  the  observed  rate  of 
transmission  loss  in  the  Arctic.  K  t  :•  ~i  •-  \ntJuJc.  !  c.r-n  * -f* c 
•  h.  '  .  K  Vi<  l  •  i  :IOvl  Jos.%. 
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Section  1 

DESCRIPTION  OF  ARCTIC  ICE  ROUGHNESS 


The  roughness  of  the  bottom  surface  of  Arctic  sea 
ice  has  been  described  using  both  discrete  and  continuous 
statistical  models.  A  discrete  model  characterizes  the 
under-ice  surface  as  a  collection  of  individual,  randomly 
distributed  ice  ridge  keels.  Probability  distributions 
describe  individual  statistics  of  the  keels  such  as  maximum 
depth,  spacing,  and  orientation.  A  discrete  model  must 
prescribe  a  representative  ridge  shape  or  an  ensemble  of 
ridge  shapes,  to  calculate  the  statistics  of  the  surface  from 
the  discrete  statistics  of  the  ice  structures.  Discrete  ice 
ridge  models  have  been  proposed  by  Diachok1,  Hibler  et 
al.2»3,  Lowry  and  Wadhams4,  Mock  et  al.5,  and  Wadhams6. 
Continuous  statistical  models  treat  the  under-ice  surface  as 
a  stochastic  process,  and  analyze  it  using  techniques  famil¬ 
iar  from  time  series  analysis.  The  surface  is  characterized 
by  its  spatial  autocorrelation  function.  Spectral  models 
give  a  more  complete  description  of  the  shape  of  the  under¬ 
ice  surface.  However,  they  are  limited  in  that  only  surfaces 
whose  depth,  h,  has  a  Gaussian  distribution  can  be  completely 
characterized  in  this  way.  Spectral  models  for  under-ice 
roughness  have  been  published  by  Ackley,  et  al.7,  Hibler 
and  LeSchack8,  Hibler9,  Rothrock  and  Thorndike,10  and 
Mellen11.  In  the  remainder  of  this  section,  probability 
distributions  for  ridge-keel  depth,  orientation,  and  spacing 
will  be  developed  from  first  principles,  spectral  models  of 
the  ice  surface  will  be  discussed,  and  techniques  for 
estimating  the  roughness  spectra  of  discrete  models  will  be 
developed. 
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1.1 


THE  DISTRIBUTION  OF  DISCRETE  ICE  STRUCTURES 


Ridge-Keel  Depth 

Using  a  simple  variational  derivation,  Hibler2 
argues  that  the  cross-sectional  area.  A,  of  the  deformed  ice 
in  an  Arctic  ice  ridge  is  a  random  variable  with  the  exponen¬ 
tial  distribution: 

1 -e_aA  .  (1) 

The  probability  density  of  the  area  is  obtained  by  differen¬ 
tiating  Eq.  (1)  with  respect  to  A.  It  is  the  familiar  expo¬ 
nential  density: 

ae"aA  .  (2) 

Hibler 's  analysis  is  based  on  a  collection  of  ridge  "links", 
ice  structures  of  fixed  length  L,  each  having  a  constant  but 
random  cross-sectional  area.  The  derivation  of  the  exponen¬ 
tial  distribution  follows  intrinsically  from  Hibler's  assump¬ 
tion  2,  that  given  a  system  of  N  ridge  links  of  length  L, 
each  with  cross-sectional  area  A,  "all  possible  ridge-link 
(area)  arrangements  yielding  the  same  net  deformation  of  a 
given  expanse  of  the  ice  sheet  are  equally  likely." 

Hibler  derives  a  distribution  for  ridge  keel  depth 
based  on  the  assumption  that  the  cross-sectional  area,  A,  is 
proportional  to  the  square  of  the  depth,  h.  This  is  a 
reasonable  assumption  if  the  shapes  of  the  ridge  cross- 
sections  are  all  similar  for  instance,  isosceles  triangles 
with  the  same  base  angle). 


A  Rayleigh  distribution  for  the  keel  depth, 

^g-abh2  ,  (3) 

results  when  the  proportionality  relation, 

A=bh 2  ,  ( 4 ) 

The  associated  probability 

(5) 

The  exponential  model  for  the  cross-sectional  area  is  con¬ 
sistent  with  a  negative  exponential  behavior  of  mean  ice 
draft  for  large  drafts,  a  behavior  observed  by  Wadhams6. 

Hibler  originally  predicted  a  modified  Gaussian 
distribution  for  keel  depth.  This  was  based  on  mistakenly 
substituting  Eq.  (4)  into  Eq.  (2),  in  the  derivation  above. 
In  a  later  publication,  Hibler3  supports  the  modified 
Gaussian  distribution  as  providing  a  good  fit  to  observed 
ridge  counts.  Lowry  and  Wadhams4  criticized  the  Rayleigh 
distribution,  when  it  was  suggested  by  Diachok1,  as  pre¬ 
dicting  too  small  a  proportion  of  small  ridges.  They  propose 
a  modification  to  the  Gaussian-type  distribution,  emphas-izing 
the  behavior  of  the  distribution  for  small  ridges.  The 
controversy  over  the  keel  depth  distribution  seems  centered 
on  the  issue  of  small  ridges. 

It  will  be  shown  in  Section  3  that  acoustic  trans¬ 
mission  loss  in  the  Arctic  depends  more  strongly  on  the 
proportion  of  large  ridges.  Because  the  asymptotic  behavior 
of  mean  ice  draft  for  large  drafts  supports  the  Rayleigh 


is  substituted  into  Eq. 
density  is: 


2abhe 


-abh2 
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model  for  large  keel  depths,  the  Rayleigh  model  was  chosen 
for  scattering-loss  calculations  in  this  report.  A  modified 
Rayleigh-type  distribution,  based  on  the  area  proportionality 
relation: 

A*bh 2  +  ch  ,  ( 6 ) 

could  be  used  to  improve  the  behavior  for  small  ridges.  The 
probability  density  arising  from  the  substitution  of  Eq.  (6) 
into  Eq.  ( 1 )  is: 

( 2abh  +  ac)  e_(abh2  +  ach)  .  (7) 
1.2  RIDGE  ORIENTATION 

It  is  commonly  assumed  for  modeling  purposes  that 
ice-ridge  orientation  is  directionally  isotropic.  This  is  a 
reasonable  assumption,  even  though  statistical  analyses  by 
Mock  et  al.5  and  Hibler9  indicate  that,  locally,  there  may  be 
a  preferred  direction  for  ridge  orientation.  However,  even 
in  the  isotropic  case,  the  collection  of  ridges  which  inter¬ 
sect  a  '’’ven  measurement  track  are  not  randomly  oriented,  but 
tend  to  lie  perpendicular  to  the  track.  The  statistics  of 
line  segments  crossing  a  given  line  are  discussed  by 
Feller12,  as  the  Buffon  needle  problem.  It  was  first 
mentioned  in  the  context  of  ice  ridge  statistics  in  Mock  et 
al.5. 

The  relative  probability  that  a  given  ridge  link  of 
length  L  will  cross  the  track  is  proportional  to  the  length 
of  the  projection  in  the  direction  perpendicular  to  the 
track,  L  sin  0,  where  9  is  the  angle  between  the  axis  of  the 
ridge  and  the  track.  The  probability  density  for  the  angle 
of  intersection  of  ridges  crossing  the  track  is: 
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( 1/2)  sin  0  d0 


r 


for  0  <  0  <  n 


9 


(10) 


with  the  associated  distribution  function: 

(1/2)  (1  -  cos  9)  .  (11) 

1.2.1  Application:  Apparent  Width  of  a  Ridge 

Any  linear  structure  of  given  width  W  will  appear 
wider,  if  it  is  crossed  at  an  angle.  The  apparent  width  is: 

W  =  W/sin  9  ,  (12) 

cl 

where  0  is  the  angle  of  intersection  between  the  axis  of  the 
structure  and  the  measurement  track.  Based  upon  the  dis¬ 
tribution  of  9  given  in  Eq.  (11),  Wa  has  the  distribution: 

/(I  -  ( W/W  )  2  on  [W,  Hr  (13) 

Q 


and  density. 


W2/(W  2  /( W  2  -  W2)  )  dw  .  (14) 

a  a  3 

An  immediate  consequence  is  that,  in  terms  of  the  actual 
width,  W,  the  mean  value  of  the  apparent  width  is: 

<W  >  ■  ( n/2)  W  .  (15) 

u 

Eq.  (14)  can  be  interpreted  as  the  conditional  probability  of 
apparent  width  for  a  given  actual  width: 

p(Wa|W)  . 
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1.2.2 


Application;  The  Width  Distribution  of  Leads  Along 
a  Track 


In  Wadhams6,  the  number  density  of  leads  per 
kilometer,  of  observed  width  d  meters,  was  found  to  obey  a 
power  law; 


n(d)  »  15  d"2  .  (16) 

A  lead  was  defined  to  be  a  continuous  sequence  of  depth 
points,  greater  than  5  meters  long  and  not  exceeding  1  meter 
of  draft,  as  observed  from  an  upward-looking  sonar. 

A  power  law  like  Eq.  (16)  is  easy  to  distinguish 
from  a  negative  exponential  law.  It  is  easy  to  devise  an 
argument,  however,  that  the  actual  width  of  leads  should  be 
exponentially  distributed.  Leads  are  formed  when  tensile 
forces  pull  the  ice  apart.  If  the  ice  opens  at  N  parallel 
leads  to  a  total  width  of  Wt  meters,  and  if  all 
combinations  of  N  lead  widths,  which  sum  to  Wt,  are  equally 
likely,  then  a  negative  exponential  density  for  lead  width 
results.  The  density  for  actual  lead  width  has  the  form: 

P(W)  *  (1/W0)  e_W/W°  .  (17) 

The  apparent  width  of  randomly  oriented  leads,  as 
observed  along  a  track,  can  be  obtained  by  "randomizing"  the 
distribution  for  actual  lead  width,  Eq.  (17),  by  the  con¬ 
ditional  distribution  for  the  apparent  lead  width,  Eq.  (14), 
in  the  form: 


This  integral  can  be  carried  out  in  closed  form,  in  terms  of 
the  special  functions  Ilf  a  Bessel  function  of  imaginary 
argument,  and  LL,  a  modified  Struve  function,  to  obtain: 

(it/(2  W0 )  )  [  I  !  '  ( Wa/W  o )  -  Lj'  (Wa/W0)]  .  (18) 

Asymptotic  expansions  for  these  functions  reveal  that  for 
large  argument,  p(Wa)  is  asymptotic  to  Wa“3.  The  mean 
apparent  width  is: 

<*/2)  W0  , 

where  W0  is  the  actual  mean  width.  Wadhams1  distribution  for 
apparent  lead  width  is  inversely  proportional  to  the  square, 
rather  than  the  cube,  of  the  width.  However,  he  speculates6 
that  in  the  region  where  the  data  was  taken,  offshore 
Greenland,  the  orientation  of  leads  and  polynyas  may  not  be 
isotropic. 

1 . 3  RIDGE  SPACING 

Ridges  appear  to  be  distributed  uniformly  along 
measurement  tracks,  see  Hibler2.  This  implies  that  ridge 
spacing  has  a  negative  exponential  distribution.  Lowry  and 
Wadhams4  have  published  a  histogram  of  observed  ridge  spacing 
which  departs  from  the  exponential  model  for  spacings  below 
about  100  meters.  Since  ridges  are  often  indistinguishable 
when  they  lie  close  together,  selection  rules  are  used  to 
decide  when  a  structure  in  the  ice  consists  of  one  ridge  or 
two.  Lowry  and  Wadhams  suggest  that  the  shortage  of  small 
ridge  spacings  in  their  histogram  is  an  artifact  of  a  selec¬ 
tion  criterion  which  counts  two  close  ridges  as  one. 
However,  to  account  for  the  entire  deviation  their  model  had 
to  use  triangular  ridges  which  are  unusually  wide  for  their 
height. 


An  alternative  explanation  for  an  actual  departure 
from  the  uniform  distribution  of  ridges  along  a  track  can  be 
based  on  the  net  expansion  of  the  ice  pack  due  to  the  ridge 
formation  process.  Ridges  form  in  cracks  in  the  ice  pack. 
The  spacing  between  these  cracks  can  be  modeled  as  a  random 
variable  with  negative  exponential  distribution.  Under 
tensile  forcing  the  cracks  open  and  are  forced  partially 
together  again,  crushing  new  ice  which  has  frozen  in  the 
opening,  to  form  the  ridge.  The  space  held  open  in  the  crack 
oy  crushed  ice  (call  it  the  "cork  width")  is  a  random  vari¬ 
able,  which  scaling  suggests  may  be  proportional  to  the 
ridge-keel  depth.  When  observed  along  a  track,  the  cork 
width  will  be  randomized  by  the  orientation  distribution. 
The  random  variable  associated  with  ridge  spacing  is,  in  this 
case,  the  sum  of  two  positive  random  variables,  the  crack 
spacing  and  the  cork  width.  The  probability  density  of  ridge 
spacing,  which  is  calculated  as  the  convolution  of  the  densi¬ 
ties  of  crack  spacing  and  cork  width,  thus  has  a  limiting 
value  of  zero  for  zero  spacing. 

Lowry's  and  Wadharas'  method  suggests  that  an 
exponential  model  of  ridge  spacing  is  reasonable,  if  the 
overlap  of  r.dges  is  taken  into  account.  The  alternative 
approach  suggests  that  the  exponential  model  has  problems  for 
small  ridge  spacing,  since  due  to  the  cork  spacing  the  ridges 
tend  to  be  a  certain  distance  apart.  For  the  purpose  of 
calculations  in  this  paper,  a  uniform  distribution  of  ridges 
will  be  assumed,  implying  an  exponential  distribution  of 
ridge  spacing.  It  will  also  be  assumed  that  mean  ridge 
spacing  is  so  large  that  the  effect  of  ridge  overlap  is 
negligible. 


1.4  SPECTRAL  PARAMETERS  OF  RIDGED  SURFACES 

The  autocorrelation  o f  a  ridged  surface  can  be 
calculated,  based  upon  the  discrete  ridge  statistics  and  a 
specified  ridge  shape  or  ensemble  of  ridge  shapes.  Coherent 
scattering-loss  formulas  and  scattering  kernels  for  the 
incoherent  field,  based  on  perturbation  techniques,  require 
the  autocorrelation  as  an  input.  The  Fourier  transform  of 
the  autocorrelation  function  is  the  roughness  spectrum  of  the 
surface.  The  autocorrelation  function  and  the  roughness 
spectrum  are  alternative  ways  of  specifying  the  second  moment 
of  the  surface.  The  specification  of  the  second  moment  is 
sufficient  information  to  completely  characterize  a  random 
surface  with  Gaussian  statistics.  A  ridged  surface  does  not 
generally  have  Gaussian  statistics  and  is  not  completely 
characterized  by  the  second  moment.  However,  a  knowledge  of 
the  second  moment  is  often  sufficient  for  predictions  of 
acoustic  scattering.  In  this  section,  techniques  for  cal¬ 
culating  the  spectrum  of  a  ridged  surface  will  be  developed, 
and  applied  to  the  example  of  randomly  oriented  cylindrical 
ridges. 

The  draft  of  the  under-ice  surface  will  be  spec¬ 
ified  by  the  depth  function, 

l  -  5(x)  , 

where  x  is  position.  The  autocorrelation  of  the  depth  is: 

W(x)  -  EU(X!  +  x)  Hxl)]/a2  ,  (19) 

where  E[  ]  indicates  the  expected  value  over  all  possible 
realizations  of  the  surface  and  a2  is  the  variance  of  the 


f" 


I 


B 


o2  -  BISIll1!  . 

It  will  be  assumed  from  this  point  that  the  surface  is 
statistically  isotropic  so  that  the  autocorrelation  depends 
only  on  the  distance  between  two  points. 

The  spectrum  is  the  Fourier  transform  of  the  auto¬ 
correlation.  There  are  one-  and  two-dimensional  spectra 
associated  with  one-  and  two-dimensional  transforms  of  the 
autocorrelation  function.  The  one-dimensional  spectrum,  Flf 
is 


F ! (K)  -  ( a2/ ( 2*) )  /  w<x)  COS  (Kx)  dx  ,  (20) 


where  K  is  the  one-dimensional  wavenumber,  in  the  form  2 it 
over  wavelength.  The  two-dimensional  transform  of  the  two- 
dimensional  isotropic  autocorrelation  function  reduces  to  a 
one-dimensional  Hankel  transform  of  the  form: 


■ 


Fj(K)  ~(o2/(2*))  /  W(x)  J0  (Kx)  K  dK  ,  (21) 

o 


where  here  K  is  the  absolute  value  of  the  two-dimensional 
wavenumber.  F^K)  is  defined  on  the  entire  real  axis, 
whereas  F2(K)  is  defined  on  the  positive  real  axis.  Occa¬ 
sionally  in  the  literature,  a  folded  version  of  Fj(K), 
defined  on  the  positive  real  axis,  is  used.  Its  value  is 
double  that  of  the  original  spectrum. 


L 


The  depth  function  £(x)  of  a  ridged  surface  is  the 
sum  of  translates  of  a  given  ridge  cross-section,  Cj(x): 


Ux) 


N 

I 

j— N 


Cj (x-Xj ) 


(22) 
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The  ridge  cross-section  itsel£  may  be  constant  or  random. 
The  variable  xj  is  the  position  o£  the  j-th  ridge. 


For  the  case  of  a  cylindrically  shaped  ridge  of 
depth  hj  and  width  R0h,  the  cross-section  is: 


/0 


Cj(X) 


Vo  , 


x  2/R  o  2 


x  <  — h  j  R  q  , 

-hjRg  <  x  <  hjR0, 
hjR0  <  x. 


(23) 


When  the  ridge  is  crossed  at  an  angle ,  it  will  appear  wider. 
The  formula  for  the  apparent  cross-section  of  the  ridge  will 
be  the  same  as  Bq.  (23),  except  that  the  width-to-height 
ratio,  Rg ,  is  replaced  by  the  apparent  ratio,  Rj .  When  the 
ridges  are  randomly  oriented,  the  apparent  ratio,  R j ,  has 
the  probability  density  given  in  Eq.  (14),  where  R0,  corres¬ 
ponding  to  W,  and  Rj,  corresponding  to  Wa,  are  the  width 
and  the  apparent  width,  respectively,  of  the  j-th  ridge, 
measured  in  units  of  height,  h. 


The  surface  with  ridges  of  elliptical  cross-section 
is  interesting  for  Arctic  applications  because  the  asymptotic 
rate  of  decay  of  its  roughness  spectrum  is  similar  to  the 
large  wavenumber  behavior  of  measured  ice-roughness  spectra. 
The  decay  rate  is  related  to  the  amount  of  small-scale  rough¬ 
ness  in  the  surface.  The  small-scale  roughness  associated 
with  an  elliptical  ridge  is  concentrated  at  the  corners  where 
the  ellipse  meets  the  flat  background.  It  is  an  open  ques¬ 
tion  whether  this  is  an  adequate  model  of  the  roughness  on  an 
ice  ridge,  or  whether  the  roughness  is  distributed  over  the 
ridge. 


The  spectrum  of  a  surface  with  ridges  of  elliptical 
cross-section  has  an  asymptotic  fall-off  of  K“3,  where  K  is 


the  wavenumber,  2x/X,  and  X  is  the  wavelength  oE  the  rough¬ 
ness  in  surface.  The  Fourier  transform  of  the  half  ellipse, 
Eq.  (22),  is  a  Bessel  function  which  has  an  asymptotic 
behavior  like  K“3/2  for  large  K.  The  asymptotic  behavior  of 
the  spectrum  is  proportional  to  the  square  of  this, 

c/K3  , 

with  some  proportionality  factor,  c. 

The  coefficient,  c,  can  be  calculated  by  Fourier 
transforming  the  autocorrelation  function,  W(x),  and  using 
the  asymptotic  properties  of  the  Bessel  functions.  Based  on 
the  definition  of  the  autocorrelation  function,  Eq.  (19),  and 
tne  stationarity  of  the  process  £,  it  follows  that: 

.  NA 

o2  W(  x )  =  lim  /  E[?(t)  £U+t)]  dt  ,  (24) 

N-*«  -NA 


where  A  is  the  mean  spacing  between  ridges.  The  spectrum  is 
found  by  substituting  Eq.  (24)  into  Eq.  (20),  where  £(x)  is 
defined  by  Eq.  (22).  The  expression  simplifies  after  inter¬ 
changing  integrals  and  Fourier  transforming  the  ridge  shapes 
to  obtain: 


Fl(K) 


E 


1  1 
Tna  7^ 


N  tih .  iKx  .  2 

l  -_1  Jl(Kh.R.)  e  3 
j  =-N  J  -> 


(25) 


where  h^ ,  R^  ,  and  x^  are  the  depth,  apparent  width  factor  and 
position  of  the  jth  ridge.  On  carrying  out  the  expected 
value  with  respect  to  ridge  position,  xj,  the  cross  terms 
in  the  absolute-value-squared  factor  in  Eq.  (25)  reduce  to 
zero,  producing: 
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F^K)  -  E  lim  ^  ji-N  T2  Jl(KhjRj)  J  •  (26) 

On  taking  the  expected  value  into  the  sum,  Eq.  (26)  simpli¬ 
fies  to: 

F  (K)  «  - —  E  [h 2  J2  (KhR)  ]  .  (27) 

24  K2  1 

The  asymptotic  behavior  of  Eq.  (27)  for  large  K  can  be  found 

in  terms  of  the  asymptotic  behavior  of  J2: 

1 

jJ(z)  ~  (1  -  sin  ( 2z)  ).  (28) 

For  large  K,  the  oscillatory  part  of  Eq.  (28)  makes  a  negli¬ 
gible  contribution  to  Eq.  (27).  The  asymptotic  behavior  of 
the  spectrum  for  large  K  reduces  to: 

F  (K)  ~  J -  E  [h]  E  { 1/R)  . 

2A  K3 

When  the  ridges  are  randomly  oriented,  the  expected  value  of 
1/R  can  be  calculated  using  Eq.  (14)  to  be:  x/(4  R0).  The 
spectrum  is  asymptotic  to: 

F»<K»  -  (ms;  E'h>)  ^  •  <29) 

The  mean  and  RMS  depth  statistics  of  the  surface 
can  be  calculated  in  closed  form.  The  mean  depth  is: 

E  [5<x>]  "  E  [2S&  /-NA  *{x)  dX^f  (30) 

where  the  stationarity  of  the  process  £(x)  has  been  used. 
Similarly,  the  RMS  depth  is: 


E  [52<x)]  -  E  52(x)  dx]. 


(31) 


When  Eqs.  (22)  and  (23)  are  substituted  for  5(x)  in  Eqs.  (30) 
and  (31 ),  the  integrals  can  be  carried  out  to  obtain: 

E  U(x)  ]  «  -jA  E  [h 2  ]  E  [R],  (32) 

and: 

E  U2(x)  ]  -  E  [h  3  ]  E  [R  ].  (33) 

For  the  case  of  randomly  oriented  ridges,  the  expected  value 
of  R  is  kRq/2 . 

In  Eq.  (33),  it  has  been  assumed  that  the  effect  of  ridge 
overlap  is  negligible. 

The  statistics  of  the  rough  surface  depend  on 
the  higher  moments  of  the  keel-depth  distribution.  For  a 
Rayleigh  distribution  of  keel  depths,  the  first  three  moments 
are: 


[h]  -  H  , 

(34) 

[h*]  -  1  K*,  and 

(35) 

[h3]  -  |  K3  . 

(36) 

Based  on  the  higher  moments  of  the  Rayleigh  keel  depth 
distribution,  and  assuming  a  random  orientation  of  ridges, 
the  derived  statistics  of  the  ridged  surface  referenced  to 
the  bottom  of  the  undeformed  level  ice,  reduce  to: 
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mean  ice  dra£t. 


EU) 


RMS  ice  draft/ 


EU2) 


(37) 


(38) 


the  standard  deviation  of  roughness  (the  variance  of  ice 
draft) , 


4R0  «Rf  2 


a -  -  —  H3  -  h4/  and 


(39) 


the  asymptotic  behavior  of  the  spectrum/ 


Fi(K>  ~1 ~ 


(40) 
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Section  2 

SAIC  INTERIM  SCATTERING  MODEL/ICE  (SISM/ICE) 


The  SAIC  Interim  Scattering  Model/Ice  (SISM/ICE)  is 
a  reflection  loss  model  for  the  coherent  reflected  field  in 
ice-covered  environments.  A  sequence  of  scattering-loss 
formulas  based  on  perturbation  techniques  provides  a  smooth 
transition  in  frequency  to  an  asymptotic  loss  level  at  the 
highest  frequencies.  The  perturbation  formulas  are  derived 
in  Bass  and  Fuks13  and  are  summarized  in  Stokes.14  The 
derivation  of  the  formula  for  the  asymptotic  loss  level  is 
contained  in  Burke  and  Twersky15  and  is  the  same  one  used  in 
Diachok.1  The  probabilistic  techniques  derived  in  Section  1 
are  used  to  relate  the  standard  deviation  of  roughness  to  the 
ice-roughness  parameters  required  by  the  individual 
reflection-loss  formulas.  The  standard  deviation  of  rough¬ 
ness  is  the  only  parameter  required  by  the  model  for  predic¬ 
tions.  A  contour  map  of  the  standard  deviation  of  roughness 
in  meters,  shown  in  Figure  1,  is  from  Buck.16  It  is  already 
available  to  Fleet  users  as  part  of  Buck's  empirical 
transmission-loss  model.  The  range  of  validity  of  SISM/ICE 
is  from  0-5000  Hz  and  0-45  degrees. 

2.1  ICE  ROUGHNESS  PARAMETERS 

The  under-ice  surface  is  considered  to  be  flat  with 
cylindrical  bosses  of  elliptical  cross  section. 

The  bosses: 

1)  have  a  half-width-to-depth  ratio  R0  of  1.6  as 
suggested  by  Diachok,1 

2)  have  a  Rayleigh  depth  distribution, 


» 


V\V 
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Figure  1.  Contours  of  standard  deviation  ice 
depth  (m)  per  LeSchack.  Dotted 
line  (...)  is  1Q00  ra  depth  contour 
and  delineates  ocean  area  for  which 
TL  models  apply. 


3)  have  random  orientation,  and 

4)  have  random  spacing  along  a  track. 

Based  on  these  assumptions,  the  standard  deviation 
of  roughness,  a,  and  the  asymptotic  behavior  of  the  spectrum 
are  related  to  the  ridge  keel  parameters  by  Eqs.  (39)  and 
(40),  where: 

R0  =»  1.6  is  the  half-width-to-depth  ratio, 

A  is  the  mean  keel  spacing,  and 

K  is  the  mean  keel  depth. 

Under-ice  roughness  spectra  measured  along  a  track 
can  be  well  approximated  by  the  two-parameter  spectral  model, 

- c - -  (41 ) 

(  p2  +  K2) 3/2 


where  K  (=2x/\)  is  the  wavenumber  of  the  roughness  component 
in  the  surface.  This  function  has  an  asymptotic  fall-off  for 
large  K  of  the  form: 


The  standard  deviation  of  roughness  derived  from  such  a  spec¬ 
trum  is: 


a2 


(43) 
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2.2 


DERIVATION  OF  ICE-ROUGHNESS  PARAMETERS  FROM  THE 
STANDARD  DEVIATION  OF  ROUGHNESS 


The  SISM/ICE  model  is  a  hybrid  of  scattering 
tneories  based  on  continuous  and  discrete  roughness  models. 
Probabilistic  techniques  are  used  to  estimate  all  the 
required  parameters  for  both  types  of  model,  based  on  the 
single  parameter,  the  standard  deviation  of  roughness.  An 
assumed  random  orientation  of  ridges  and  effective  ridge 
shape  are  used  to  calculate  a  surface  roughness  spectrum. 
Two  discrete  parameters,  the  mean  ridge  spacing  and  depth, 
are  calculated  as  inputs  to  the  Burke  and  Twersky15 
scattering  theory. 

If  the  mean  ridge  spacing  is  known,  the  mean  keel 
depth  can  be  estimated  from  the  standard  deviation  of  rough¬ 
ness  by  solving  Eq.  (39)  as  the  limit  of  the  following 
iterative  algorithm: 

hn+1  =  [(A/4R0)(o2  +  it2  R2  h£/A2)  ]1/3.  (44) 

For  the  purposes  of  this  model,  it  is  assumed  that  the  mean 
ridge  spacing  is: 

A  =  100  m. 


An  improved  estimate  of  A  on  a  regional  basis  could  be  used 
to  improve  estimates. 


The  spectral  scattering  models,  based  on  the  spec¬ 
tral  model,  Eq.  (41),  require  the  two  parameters. 


c 


uh 

8  AR  o  ' 


(45) 
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from  Eq.  (40),  and: 


P  *  ( 2c/ a2) 1/2  ,  (46) 

from  Eq.  ( 43 ) . 

The  autocorrelation  function  related  to  Eq.  (41)  is: 
pr  Kj  (pr), 

where  KL  is  a  Bessel  function  of  imaginary  argument  and  r  is 
the  separation  distance.  The  two-dimensional  isotropic 
roughness  spectrum  related  to  Eq.  (41)  is: 

2c 

( P2  +  K2) 2 

2.3  REFLECTION  LOSS  FORMULAS 

In  terms  of  the  ice-roughness  parameters,  four 
formulas  for  the  reflected  power  |  R2 1  are  evaluated.  The 
predicted  value  of  | R2 |  is  the  maximum  of  the  four  calculated 
values. 

•  Low  Frequency-Free  Surface  Formula 

| R |  2  *  1  -  4  (sin9)  k2a2(  1-it2c  (sin0)/2),  (47a) 

•  High  Frequency-Free  Surface  Formula 

|  R 1  2  »  1  -  4  (sine)  (1.198)  c  ( 1c/  p )  3/2  ,  (47b) 

•  High  Frequency-Rigid  Surface  Formula 
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x  <  sine 


9 


I  R| 2 


(sin8  -  x)  2  +  x2 
(sine  +  x)2  +  x2 
.2 


x  >  sine  , 


(47c) 


where: 

x  =  1.311  c  ( k/ 0 ) 1  /  2  . 

•  Asymptotic  Twersky  formula 

| R | 2  *  (sine  -  x)2/(sine  +  x)2  ,  (47d) 

where : 

x  n  L  cos(  6-y)  , 
tany  *  p“2tan0, 

L2  =•  w2( tan2  e  +  p4)/( tan2  e  +  p2), 
w  »  Rg  n  h/2, 

P  *  2/(RqTc), 

n  ■  1/A. 

In  all  of  formulas  47a-47d,  6  is  grazing  angle  and  k  is  the 
acoustic  wavenumber  in  inverse  meters. 

2.4  DISCUSSION 

SISM/ICE  is  suggested  as  an  interim  scattering-loss 
model  in  ice-covered  environments.  It  is  based  on  several 
simplifications  and  assumptions.  However ,  it  fulfills  the 
requirement  for  a  loss  model  with  a  wide  range  of  validity, 
based  on  the  only  known  environmental  parameter,  the  standard 
deviation  of  roughness.  SISM/ICE  can  be  used  with  acoustic 
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propagation  models  such  as  FACT  and  RAYMODE  to  predict  the 
coherent  propagating  field.  However,  ray-based  models  may 
not  adequately  reproduce  source  and  receiver  depth  effects 
due  to  modal  propagation  at  the  lower  frequencies.  There  may 
also  be  a  significant  incoherent  component  propagating  at  the 
higner  frequencies,  which  would  not  be  predicted  by  these 
models.  As  propagation  models  become  available  which  can 
describe  the  propagating  incoherent  field,  SISM/ICE  can  be 
extended  to  give  the  scattering  kernel  for  the  incoherent 
field.  The  details  for  doing  this  are  contained  in  Refs.  13, 
14,  and  15. 


2-7 


Section  3 

HIGH-ANGLE  PE  SIMULATION  OF  ARCTIC  TRANSMISSION  LOSS 
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The  High-Angle  Parabolic  Equation  (HAPE)  has  been 
used  to  simulate  transmission-loss  data  reported  by  Diachok1, 
taken  at  two  locations  in  the  Beaufort  Sea.  The  data  sets 
are  representative  of  the  Central  Arctic  environment  and  the 
off-shore  Canadian  environment.  The  dominant  loss  mechanism 
was  found  to  be  scattering  from  Arctic  sea-ice  ridge  keels. 
For  the  purpose  of  the  HAPE  simulations,  equally  spaced 
trapezoidal  ridges  were  used.  Two  simulated  transmission 
loss  runs  were  carried  out  at  each  location,  the  first  with 
random  depth  keels,  the  second  with  equal  sized  keels  whose 
depth  was  equal  to  the  mean  depth  of  the  random  keels.  For 
both  data  sets,  the  random  depth  keels  produced  significantly 
more  loss  than  the  constant  depth  keels,  and  gave  good  agree¬ 
ment  with  the  data. 

3.1  ENVIRONMENTAL  DATA 

The  random  keel  depths  used  in  the  simulations 
satisfy  the  Rayleigh  distribution,  Eq.  (5).  This  is  a  one- 
parameter  model,  dependent  on  the  mean  keel  depth.  In  the 
Central  Arctic  case,  the  mean  keel  depth  and  the  mean  ridge 
spacing  suggested  by  Diachok1  were  used  in  the  simulation. 
In  the  off-shore  case,  the  mean  keel  depth  was  estimated  from 
a  histogram  of  ridge-keel  depths  published  by  Hibler2.  These 
data  are  based  on  upward-looking  sonar  measurements  carried 
out  by  the  USS  Sargo  in  the  winter  of  1960,  near  the  location 
of  the  off-shore  transmission-loss  experiment.  Hibler's  data 
provide  statistics  for  keel  depth  and  mean  spacing  of  ridges 
deeper  than  6.1  meters.  The  keel-depth  distribution  for 
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smaller  ridges  has  been  extrapolated,  based  on  a  least 
squares  £it  of  the  Rayleigh  distribution  to  the  histogram. 
In  addition,  the  overall  mean  ridge  spacing  has  been  esti¬ 
mated  from  the  published  figure,  which  once  again  relates  to 
ridges  deeper  than  6.1  meters,  by  multiplying  by  the  esti¬ 
mated  percentage  of  ridges  deeper  than  6.1  meters.  The  width 
of  the  ridges  was  assumed  to  be  3.2  times  the  height,  consis¬ 
tent  with  Diachok.1  In  the  HAPE  simulations,  the  ridge  width 
was  increased  by  a  factor  of  n/2,  to  account  for  the  mean 
apparent  expansion  due  to  random  orientation. 

In  all  simulation  runs,  the  compressional  sound 
speed  in  the  ice  was  assigned  the  value  2280  m/s.  The 
density  was  set  at  .89  g/cc.  The  sound  speed  profile  in  the 
water  was  linearly  interpolated  in  depth  between  the  values 
in  the  following  table: 


depth  (m) 

speed  (m/s) 

0.0 

1436.0 

50.0 

1436.0 

100.0 

1443.0 

300.0 

1457.8 

700.0 

1466.6 

1085.0 

1473.5 

A  water  depth  of  1085  meters  was  used  in  all  simu¬ 
lations.  The  ocean  bottom  was  assumed  to  be  totally 
absorbing.  In  the  off-shore  case,  the  measurement  tracks  lie 
over  shallow  ridge  and  abyssal-hill  bottom  environments, 
associated  with  the  Lomonosov  ridge  and  Alpha  Cordillera. 
The  bottom  loss  for  the  off-shore  case  can  be  expected  to  be 
quite  high.  In  the  Central  Arctic  case,  the  measurement 
tracks  lie  over  the  Canada  Abyssal  Plain,  where  water  depth 
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averages  approximately  3500  meters.  The  ocean  bottom  in  the 
abyssal  plain  environment  can  be  expected  to  be  a  good 
reflector.  Buck18  observes  that  in  the  Beaufort  Sea,  the 
bottom-reflected  energy  makes  up  the  major  part  of  the  total 
received  energy  for  ranges  up  to  55  km. 

3.2  ANALYSIS  OF  DATA  AND  CONCLUSIONS 

The  comparison  between  HAPE  simulations  and  data 
from  the  Central  Arctic  environment  is  shown  in  Figures  2  and 
3.  Trapezoidal  ice  ridges  with  a  mean  depth  of  4.3  meters 
and  spacing  of  9.5  ridges  per  kilometer  were  used.  These  are 
the  values  suggested  by  Diachok1,  in  his  analysis  of  the 
data.  In  Figure  2,  the  ridges  have  a  constant  depth  and  are 
equally  spaced.  At  both  40  and  50  Hz,  the  simulation  under¬ 
estimates  the  observed  loss  at  all  ranges.  At  300  km,  the 
difference  is  about  5  dB.  In  Figure  3,  the  ridges  have  a 
random  depth  with  Rayleigh  distribution.  The  simulation  with 
random  depth  ridges  gives  good  agreement  with  the  data  at 
long  ranges,  while  overestimating  the  loss  at  short  ranges. 
Buck18  demonstrates  convincingly  that  the  differences  at  the 
short  ranges  are  due  to  bottom-reflected  energy,  which  was 
not  included  in  the  simulation.  The  implications  of  this 
numerical  experiment  for  future  modeling  efforts  are: 

1 )  Scattering  from  random  depth  ice-ridge  keels 
can  account  for  observed  transmission-loss 
levels  in  the  Arctic. 

2)  Random  depth  keels  with  a  Rayleigh  depth  dis¬ 
tribution  produce  more  transmission  loss  than 
equal  depth  keels  of  the  same  mean  depth. 

The  relevant  ridge  statistics,  mean  keel  depth 
and  spacing,  required  for  transmission-loss 
calculations  can  be  estimated  from  observa¬ 
tions. 
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Figure  2 


Based  on  a  simulated  ice  surface  con¬ 
sisting  of  equal  depth,  equally  spaced 
trapezoidal  ice  keels,  ransmission-loss 
data  at  the  Central  Arctic  site  at  40  Hz 
(....)  and  50  Hz  (xxxx)  is  compared  with 
simulated  transmission  loss  at  40  Hz 
(  )  and  50  Hz  {----),  using  the 
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Figure  3.  Based  on  a  simulated  ice  surface  con¬ 
sisting  of  random  depth,  equally  spaced 
trapezoidal  ice  keels,  transmission-loss 
data  at  the  Central  Arctic  site  at  40  Hz 
(....)  and  50  Hz  (xxxx)  is  compared  with 
simulated  transmission  loss  at  40  Hz 

( - )  and  50  Hz  ( - ),  using  the 

Arctic  High  Angle  PE. 


4)  Scattering  is  dominated  by  the  deeper  ridges, 
implying  that  ridge  observations  should  empha¬ 
size  the  larger  ridges. 

5)  The  bottom-bounce  path  may  produce  an  important 
contribution  at  short  ranges  in  an  abyssal 
plain  environment. 

The  transmission-loss  data  from  the  measurement 
tracks  in  the  off-shore  Canadian  environment  are  shown  in 
Figure  4  at  40  Hz  and  Figure  5  at  50  Hz,  together  with  HAPE 
simulations.  A  mean  keel  depth  of  7.7  meters  and  keel 
spacing  of  10  per  kilometer  were  estimated  from  Hibler's2 
histogram  of  observed  keel  depths  for  the  Offshore  Province, 
subprovince  1 ,  located  near  the  middle  of  the  measurement 
tracks.  These  keel  statistics  correspond  to  a  standard 
deviation  of  roughness  of  about  4.5  meters. 

In  Figures  4  and  5,  the  solid  line  corresponds  to 
the  simulation  results  with  equal  sized  ridges,  the  dashed 
line  to  random  depth  ridges  of  Rayleigh  depth  distribution. 
Once  again,  the  random  depth  keels  produce  more  loss  than  the 
equal  depth  keels,  in  this  case  on  the  order  of  10  dB  at  a 
range  of  300  km.  The  random  depth  simulation  mimics  the 
slope  of  the  data  extremely  well,  but  underestimates  observed 
loss  by  a  constant  4  dB  at  40  Hz  and  7  dB  at  50  Hz.  The 
40-Hz  data  was  shifted  by  3  dB  in  an  unknown  direction  as 
published  by  Diachok,  suggesting  that  the  difference  may  be 
7  dB  in  both  cases. 

Despite  the  differences  in  level,  several  con¬ 
clusions  are  possible: 

1)  Random  depth  keels  with  a  Rayleigh  depth 
distribution  produce  more  transmission  loss 
than  equal  depth  keels  of  the  same  mean  depth. 
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Figure  4.  Simulated  transmission  loss  at  40  Hz, 
using  the  Arctic  High  Angle  PE,  based  on 
a  simulated  off-shore  Canadian  ice  ridge 
environment  with  equal  depth,  equally 

spaced  trapezoidal  keels  ( - )  and 

random  depth,  equally  spaced  trapezoidal 

keels  ( - )  are  compared  with  40  Hz 

transmission-loss  data. 
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Figure  5.  Simulated  transmission  loss  at  50  Hz, 
using  the  Arctic  High  Angle  PE,  based  on 
a  simulated  off-shore  Canadian  ice  ridge 
environment  with  equal  depth,  equally 

spaced  trapezoidal  keels  ( - )  and 

random  depth,  equally  spaced  trapezoidal 

keels  ( - )  are  compared  with  50  Hz 

transmission-loss  data. 
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2)  Scattering  from  random  depth  ice-ridge  keels 
with  Rayleigh  depth  distribution  can  account 
for  the  observed  rate  of  decay  in  range  of 
observed  transmission-loss  levels  in  the 
Arctic. 

3)  The  ridge  statistics,  mean  keel  depth  and 
spacing,  supporting  the  transmission— loss 
calculations  can  be  estimated  from  observa- 
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